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Abstract. The current state-of-the-art of population synthesis is reviewed. The field is currently 
undergoing major revisions with the recognition of several key processes as new critical ingredi- 
ents. Stochastic effects can artificially enhance or suppress certain evolutionary phases and/or 
stellar mass regimes and introduce systematic biases in, e.g., the determination of the stellar 
initial mass function. Post-main-sequence evolution is often associated with irregular variations 
of stellar properties on ultra-short time-scales. Examples are asymptotic giant branch stars and 
luminous blue variables, both of which are poorly treated in the models. Stars rarely form in 
isolation, and the fraction of truly single stars may be very small. Therefore, stellar multiplicity 
must be accounted for since many systems will develop tidal interaction over the course of their 
evolution. Last but not least, stellar rotation can drastically increase stellar temperatures and 
luminosities, which in turn leads to revised mass-to-light ratios in population synthesis models. 
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1. Prom Fundamental Data to Cosmology 

Stellar population synthesis is at the heart of spectral energy distribution (SED) studies 
of galaxies. The very nature of population synthesis makes it a rather multi-disciplinary 
research area. It is built on our knowledge of stellar astrophysics, which in turn relies on 
fundamental data derived from laboratory physics, such as atomic line data and nuclear 
reaction rates. On the other hand, those studying galaxy evolution trust population 
synthesis models for the interpretation of galaxy colors and their evolution with cosmic 
time. Ultimately, cosmological models often rest on our faith in population synthesis. 
The chain of assumptions that must be made is usually — subconsciously — assigned 
errors which hierarchically decrease from the fields of cosmology over stellar astrophysics 
to laboratory data. One goal of this review is to raise the awareness of radically new 
developments and uncertainties in hitherto considered "well-known" areas and to foster 
feedback between the stellar and cosmological communities. 

2. The Basics of Population Synthesis 

Population synthesis aims to reproduce the observed galaxy SED from a set of pre- 
specified input parameters. Evolutionary spectral synthesis is a special, widely used case 
of population synthesis, whose goal is to reproduce observed spectra self-consistently from 
the star-formation history of a galaxy and from stellar evolution models. This method 
was pioneered by Tinsley (1980) after the advent of modern stellar evolutionary tracks 
in the Hertzsprung- Russell diagram (HRD). Evolutionary synthesis has relatively few 
free parameters but its success depends on the reliability of the adopted stellar models. 
Therefore, assessing these assumptions becomes critical. 
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The astrophysical ingredients entering evolutionary synthesis are as follows: (i) Quan- 
tities related to the star- formation process, i.e., the star- formation rate and its evolution 
with time and the stellar mass spcctriim at birth, also known as the initial mass function 
(IMF). See, e.g., the conference proceedings of Treyer ct al. (2011). (ii) Once stars have 
formed, stellar evolution models provide a prescription for the variation of the stellar 
luminosity (L), effective temperature (Tcff), and mass (M) as a function of chemical 
composition, initial mass, and time, (iii) Spectral libraries describe the spectrum of each 
star for any given {L, T^s, M). (iv) Second-order effects, such as dust attenuation or 
geometric effects are often accounted for in a very approximate way or even neglected 
altogether. 

Several widely used synthesis models are made available to the community via ded- 
icated web sites. Among the most popular packages are GALEV (Schulz et al. 2002; 
Kotulla et al. 2009), GALEXEV (Bruzual & Chariot 1993; 2003), PEGASE (Fioc & 
Rocca-Volmerange 1997; Le Borgne et al. 2004), and Starburst99 (Leitherer et al. 1999; 
Leitherer & Chen 2009). In general, the predictions made by different synthesis codes 
are in good agreement. When discrepancies are found, they can usually be understood 
in terms different intrinsic input parameters and/or treatment of peculiar stellar evolu- 
tionary phases, such as the asymptotic giant branch (AGB) phase (Vazquez & Leitherer 
2005; Conroy & Gunn 2010). 

The mentioned synthesis packages restrict their output to models for the stellar lumi- 
nosities, with some generic contribution by nebular emission. More specialized applica- 
tions generate panchromatic galaxy SEDs accounting for the stellar as well as the gas 
emission using full photo-ionization modeling (Dopita et al. 2005; 2006a; 2006b; Groves 
et al. 2008). SEDs with self-consistent inclusion of dust were first introduced by Dwek & 
Scalo (1980) and refined most recently by Dwek & Cherchneff (2011). Chemical evolution 
models for the full life-cycle of the major elements over the cosmic evolution of galaxies 
were discussed by Mattoucci (2009; 2010). 

Rather than reviewing these "traditional" synthesis models I will focus on some re- 
cent, new developments which may challenge well-established results: stochasticity, pe- 
culiar evolutionary phases, stellar multiplicity, and stellar rotation. The importance of 
these mechanisms and concepts had been recognized for quite some time but only in 
recent years have they been implemented in population synthesis codes for quantitative 
comparison with observations. 

3. Stochasticity 

Traditional population synthesis models scale all population properties by the star- 
formation rate. This approach is strictly valid only in the limit of a stellar population 
containing an infinite number of stars. Real star clusters or galaxies, however, contain a 
finite number of stars. Furthermore, most of the light is sometimes provided by very few 
bright stars, in particular in the near-infrared (IR). The so-called stochastic fluctuations 
in the integrated population properties are the result of the random presence of these 
luminous stars. Depending on population age and wavelength domain, stochastic effects 
can become important for masses as high as ~105 Mq (Bruzual 2002; Cervino et al. 
2002; Fouesneau & Langon 2010). 

da Silva et al. (2011) developed and released SLUG, a new code to "Stochastically 
Light Up Galaxies" . SLUG synthesizes stellar populations with a Monte Carlo technique 
to treat stochastic sampling properly. Included are the effects of clustering, the stellar 
IMF, star-formation history, stellar evolution, and cluster disruption. Their code has 
been built using the same ingredients as in Starburst99, thereby allowing the immediate 
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Figure 1. R-band, far-ultraviolet (UV), and ionizing luminosities vs. time for galaxies with 
constant star-formation rates of 1 (left) and 0.01 (right) Mq yr~^ assuming clustered star 
formation. Starburst99 models (solid black lines) are compared with simulations from SLUG. 
The SLUG models are represented by their mean (dash-dotted line), median (dashed line) and 
5-95 percentile range (shaded region). From da Silva et al. (2011). 



separation of stochastic effects. An example is shown in Fig. [TJ This figure highlights 
the increased scatter and systematic offsets due to stochasticity and clustering at low 
star-formation rates. 

GALEX UV imagery suggests that star-formation rates inferred from Ha in galactic 
environments characterized by low stellar and gas densities tend to be less than those 
based on the UV luminosity (Lee et al. 2009). The origin of the discrepancy is still under 
debate. One possible explanation is that the stellar IMF is systematically deficient in 
high-mass stars in such environments. However, effects of uncertainties in the stellar 
evolutionary tracks and model atmospheres, uncertain metallicities, non-constant star- 
formation histories, leakage of ionizing photons, departures from Case B recombination, 
dust attenuation, and finally stochasticity in the formation of high-mass stars need to be 
considered. 

Eldridge (2011) found that the scatter and variation of the Ha- and UV derived star- 
formation rates are less dependent on the IMF but more on the star-formation history 
of each individual galaxy. The general trend is that those systems with weaker star- 
formation in the last 10 Myr have relatively lower Ha rates, while those with most of 
the star formation in the last 10 Myr display the opposite behavior even at low mean 
star-formation rates. Therefore any simulation attempting to predict the properties of a 
sample of galaxy must take into account the stochastic nature of star-formation. If there 
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Figure 2. Comparison of observed and modeled near-IR colors as a function of age. The data 
are derived from star clusters in the Magellanic Clouds. The models include the predictions from 
Maraston (2005; dot-dashed lines), Bruzual & Chariot (2003; triple-dot-dashed lines) and two 
FSPS models (dashed and solid lines). The observed data are individual measurements (squares) 
and averages over many star clusters (circles). From Conroy & Gunn (2010). 

are only a few clusters, the appearance of the galaxy-wide stellar population may be very 
different from what might be expected for a simple stellar population with a smooth 
star-formation history. 



4. Rapid Evolutionary Phases 

Rapid evolutionary phases pose particular challenges for population synthesis. Ther- 
mally pulsing AGB stars are especially notorious (Maraston 2011). Their uncertainties in 
stellar evolution modeling arise from short pulsation time-scales (~10* yr), the double- 
shell burning, and strong mass loss in this phase. The spectral modeling during the 
carbon-rich phase adds to the complexity. Yet AGB stars can often make a significant 
contribution to the galaxy SED in the optical and near-IR because the AGB phase is 
bright and occurs around '^1 Gyr when the overall stellar population has already faded. 
Different population synthesis models disagree quite severely in their predicted lumi- 
nosity contribution from thermally pulsing AGB stars. The Flexible Stellar Population 
Synthesis (FSPS) technique of Conroy et al. (2009; 2010) and Conroy & Gunn (2010) 
is particularly well suited for addressing this issue. FSPS is a novel Monte-Carlo-type 
approach that is capable of flexibly handling various uncertain aspects of stellar evolu- 
tion such as, e.g., thermally pulsing AGB stars, the horizontal branch, or blue stragglers. 
The technique allows one to marginalize over uncertain evolutionary aspects in order 
to understand the full uncertainties associated with the physical properties of galaxies 
including stellar masses, ages, metallicities, and star-formation rates. An example is in 
Fig. [5] where the AGB phase as implemented in different evolutionary tracks is tested. 
The conclusion is that models generally overestimate this phase, which is consistent with 
the result of Kriek et al. (2010) who used a sample of ^60 post-starburst galaxies at 
z « 1.6 in the NEWFIRM medium-band survey to assess different models. The AGB 
contribution to the integrated SED is a factor of ^3 lower than predicted. This could be 
due to lower bolometric luminosities, shorter lifetimes, and/or heavy dust obscuration of 
AGB stars. 

AGB stars affect population synthesis models both directly (via their luminosity contri- 
bution) and indirectly (via their strong mass loss influencing stellar evolution) . Luminous 
blue variables (LBVs) mirror AGB stars in their eruptions and their high mass-loss rates. 
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Figure 3. H/3 equivalent widtli vs. time for a series of instantaneous burst models generated 
with the BPASS code. Solid and dashed lines indicate populations without and with close binary 
stars, respectively. Both model families with and without binaries are shown for five different 
metallicities, higher metallicities generally having higher equivalent widths. From Eldridge & 
Stanway (2009). 

They introduce uncertainties in population synthesis similar to those of AGB stars. They 
generally do not contribute significantly to the galaxy SED but their mass ejections are 
decisive for the final evolution of massive stars. It has recently been recognized that most 
mass loss in massive stars prior to the Wolf-Rayet phase occurs in LBVs, whereas steady 
stellar winds close to the main-sequence are less important (Smith 2010). LBVs belong to 
a diverse class of transients in the upper HRD spanning a wide range of outburst luminos- 
ity which ranges from that of classical novae to that of core-collapse supernovae (Smith 
et al. 2011). The ubiquity of eruptive phases and highly time- variable evolutionary phase 
introduce formidable challenges to population synthesis. 

5. Stellar Multiplicity 

Massive stars are known for their high-degree of multiplicity. Most OB stars are found 
in binaries and multiple systems. Combining information from these various surveys, 
Mason et al. (2009) found a minimum multiplicity fraction close to 70%. Even bona fide 
single field stars may have been part of a multiple system in the past, then ejected by 
a supernova explosion or by dynamical interaction. While multiplicity deserves studies 
by its own right, neglecting this fundamental property can be a serious omission in 
population synthesis (Vanbeveren 2010). In the following we will not focus on the well- 
known vertical main-sequence displacement in the color-magnitude diagram caused by 
binaries but rather on the modified stellar evolution of close binary system and their 
consequences for population synthesis. Relevant in this context is the fraction of close 
binaries, i.e., those binaries whose evolution differs from that of single stars. Sana & 
Evans (2011) estimate a spectroscopic binary fraction of close to 50%, with a significant 
fraction likely to be interacting. 

de Mink et al. (2009) demonstrated that the evolution of massive, close binaries still 
on the main-sequence sharply differs from single-star evolution. Tidal effects can lead to 
significant spin-up and high rotational velocities. Mixing processes induced by rotation 
may be so efficient that helium produced in the center is mixed throughout the envelope. 
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Such stars evolve almost chemically homogeneously. If the metallicity is low, they remain 
blue and compact, while they gradually evolve into Wolf-Rayet stars and possibly into 
progenitors of long 7-ray bursts. 

Alternatively, at low rotation velocities, the distinguishing property of binary evolution 
will be Roche-lobe overflow, which typically occurs off the main-sequence. Depending 
on the system properties, the primary will transfer material to the secondary, which 
will accrete material and gain mass. The crucial point for population synthesis is the 
"rejuvenation effect" (van Bever & Vanbeveren 1998; Eldridge & Stanway 2009): the 
previously less massive secondary will become more massive and more luminous and will 
therefore appear younger. This is illustrated in Fig. |3l which shows the evolution of the 
nebular H/3 emission-line equivalent width of single stellar populations. The models were 
obtained with the BPASS (Binary Population and Spectral Synthesis) code of Eldridge et 
al. (2008), which accounts for the evolutionary effects of massive binaries. Compared with 
single-star models, binary populations are predicted to have larger H/3 equivalent width 
at older ages due to the rejuvenation effect. H/? is commonly used as an age indicator 
in H II regions (Stasiriska et al. 2001) and would require recalibration. For instance, 
in single-star models, a large H/3 value at older ages would be incompatible with an 
instantaneous burst and suggest ongoing star formation. Binary models for single stellar 
populations, however, might still produce large enough H/3 emission at older ages. 

6. Stellar Rotation 

Massive stars exhibit significant helium and nitrogen enrichment on the surface while 
still on the main-sequence (Hunter et al. 2009). These telltale signs of nuclear processing 
require an efficient transportation mechanism to expose materials from the convective 
core. While mass removal by stellar winds had previously been thought to assume this 
role, the main-sequence mass-loss rates are now considered too low (see Section 4) to make 
this a viable option. Alternatively, mixing, mass transfer, and mass loss in close binaries 
might cause the enrichment. However, it seems unlikely that the evolution of all massive 
stars is driven by binarity. Many stars with surface enrichment are also fast rotators. 
While the relation between nitrogen overabundance and rotation velocity is complex and 
far from unique, there is consensus that rotation is an important, if not the prime driver 
of mixing in massive stars (Brott et al. 2011a). Stellar evolution with rotation has been 
modeled over the past decade, but it is only now that quantitative evolutionary tracks 
with rotation have become available for implementation of population synthesis (Brott 
et al. 2011a; Ekstrom et al. 2011). 

Fig.|4]illustrates the principal effects caused by rotation. At high masses (M > 2OM0), 
the dominant mechanism is the increase of the convective core and the reduced surface 
opacity due to the He enhancement, resulting in a higher L and higher Teg, respectively. 
(Recall that H is the major opacity source and any decrease of its relative abundance by 
mixing lowers the opacity and therefore increases Teff.) Therefore tracks with rotation 
tend to fall above and to the left of non-rotating tracks. At lower masses (e.g., the track 
for M — 15 Mq in Fig. 2]), these effects become less important and are surpassed by 
the impact of the centrifugally induced radius increase and the associated lower Teff. At 
masses below ^2 Mq no rotational effects are expected since low-mass stars arrive on 
the zero-age main-sequence with essentially no velocity because of magnetic braking. 

We implemented the new suite of 48 different stellar evolutionary tracks, both rotating 
and non-rotating, by Ekstrom et al. (2011) in Starburst99. The grid has solar chemical 
composition and covers the full mass range from 0.8 to 120 Mq. The rotating models 
start on the zero-age main-sequence with a rotation velocity of 40% of the equatorial 
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Figure 4. Evolutionary tracks for rotating stars with initial masses of 15, 25, and 50 M© 
having solar chemical composition. Line types indicate different initial rotation velocities on the 
zero-age main-sequence. From Brott et al. (2011b). 



break-up velocity. Guidance on the choice of the initial rotation velocity is provided by 
measurements in B stars but data for O stars are still scarce. Clearly this parameter 
must be considered with care. The evolution is computed until the end of the central 
carbon-burning phase, the early AGB phase, or the core helium-flash for the massive, 
intermediate, and low-mass stars, respectively. Vazquez et al. (2007) used an earlier, 
preliminary version of these models to test the impact of stellar rotation on synthetic 
population properties. The new model suite will allow us to verify and extend the previous 
work (Levesque et al. 2012, in preparation). 

Rotation increases both L and Tcff for stars more massive than '--^20 Mq. This effect 
sets in very early on the main-sequence. As a result, a stellar population containing 
hot, massive stars is predicted to be more luminous for a given mass, and its SED is 
shifted to shorter wavelengths. Therefore the most dramatic changes with respect to 
prior models occur at the short-wavelength end of the SED. Examples are reproduced 
in Fig. [SJ which compares Mbo\ and the luminosity of the nebular Br7) resulting from 
models with and without rotation. The single stellar population models with rotation for 
-^^Boi become more luminous by ~0.4 mag (left part of Fig. [5]). The right part of Fig. [5] 
shows the Br7 luminosity for an equilibrium population. Br7 is often used to determine 
star-formation rates in obscured galaxies, such as ultra-luminous IR galaxies. For this 
specific example, the rates following from the models with and without rotation would 
be 100 and 175 M© yr~^, respectively. The extreme UV increases by a factor of a few in 
the hydrogen ionizing continuum and by several orders of magnitude in the neutral and 
ionized helium continua with the new models. The predictions for the latter need careful 
testing, as the photon escape fraction crucially depends on the interplay between the 
stellar parameters supplied by the evolution models and the radiation-hydrodynamics 
of the atmospheres. In contrast, the escape of the hydrogen ionizing photons has little 
dependence on the particulars of the atmospheres and consequently is a relatively safe 
prediction. 
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Figure 5. Comparison of stellar population properties resulting from evolutionary tracks with 
and without rotation. Left: bolometric luminosity vs. time for a single stellar population with a 
mass of 10^ Mq. Right: Br7 luminosity for continuous star formation at a rate of 100 Mq yr~^. 
Solid lines: new tracks with rotation; dashed: new tracks with zero rotation velocity; dotted: 
previous generation of Geneva tracks released 1992 - 1994. The models are for solar chemical 
composition and use a standard Kroupa IMF between 0.1 and 100 Mq. 

7. The New Standard 

Stochasticity, exotic stellar evolutionary phases, multiplicity, and rotation are now 
considered essential ingredients in population synthesis models. A new suite of population 
synthesis models which account for these effects and provide quantitative predictions 
have become available in the very recent past. Somewhat frustratingly, implementation of 
these processes cannot rely on first principles, which essentially adds new free parameters 
requiring empirical calibration. Care is therefore required when the first model generation 
is used to analyze galaxy SEDs. The community is encouraged to perform extensive 
testing and help making these new synthesis models the new standard. 
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